Late Time Neutrino Masses, the LSND Experiment and the Cosmic Microwave 

Background 
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Models with low-scale breaking of global symmetries in the neutrino sector provide an alterna- 
tive to the seesaw mechanism for understanding why neutrinos are light. Such models can easily 
incorporate light sterile neutrinos required by the LSND experiment. Furthermore, the constraints 
on the sterile neutrino properties from nucleosynthesis and large scale structure can be removed 
due to the non-conventional cosmological evolution of neutrino masses and densities. We present 
explicit, fully realistic supersymmetric models, and discuss the characteristic signatures predicted 
in the angular distributions of the cosmic microwave background. 



Introduction — The LSND experiment found evidence 
for the oscillations D„ — > v„ and v., — ► v e with an oscilla- 



tion probability of around 3x 10 



and a Am 2 > leV 



The statistical evidence for the anti-neutrino oscillations 
is much stronger than for the neutrino case, with some 
analyses finding a 5<r effect |2( . While other experiments 
restrict the regions of parameter space that could explain 
the LSND data, they do not exclude the LSND result Q . 

The confirmation of solar and atmospheric neutrino os- 
cillations has led to a "standard" framework for neutrino 
masses, with three light neutrinos resulting from the see- 
saw mechanism and the heaviest right-handed neutrino 
not far from the scale of gauge coupling unification [^|. 
The LSND result conflicts with this framework, and, 
if confirmed by the Mini-BooNE experiment [5|, would 
throw neutrino physics into a revolutionary phase. There 
are three major challenges to incorporate the LSND re- 
sult into the standard framework for neutrino masses. 
A fourth light neutrino is needed with mass in the eV 
range, and, given the Z width, this neutrino must be 
sterile. Such states are anathema to the seesaw mecha- 
nism. Secondly, neutrino oscillations in the early universe 
would ensure that this fourth neutrino state is thermally 
populated during the big bang nucleosynthesis (BBN), 
changing the expansion rate and yielding light element 
abundances in disagreement with observations 0. In 
fact, theories with four neutrino species give very poor 
fits to the combined data of the LSND experiment and 
oscillation experiments with null results, while fits in the- 
ories with 5 neutrinos are much better Q. This would 
imply that A^rbn > 5, in gross disagreement with ob- 
servations H, [lj ■ Finally, the combination of large scale 
structure surveys and WMAP data [ToL Hl| has led to a 
limit on the sum of the neutrino masses of about 0.7 eV, 
which is significantly less than the best fit values for the 
LSND neutrino masses 0, El- 

An alternative explanation for why the neutrinos are 
light has been explored recently: the scale of neutrino 
masses can be dictated by a low scale / of breaking 
of global symmetries in the neutrino sector 



characteristic feature of this mechanism is that, cosmo- 
logically, neutrinos remain massless until the symmetry- 
breaking phase transition, which occurs quite late in the 
history of the universe — hence the name "late time neu- 
trino masses". In this letter, we argue that, unlike the 
traditional seesaw framework, this alternative scenario 
can easily accommodate the sterile neutrinos required by 
the LSND experiment. Moreover, the cosmological evo- 
lution of neutrino masses and densities in this scenario 
is non-standard and, as a result, the apparent contra- 
diction between the parameters preferred by LSND and 
cosmology can be avoided, provided that / is of order 
100 keV. (We will show that such values of / can arise 
naturally in supersymmetric theories.) At the same time, 
the scenario predicts potentially observable characteristic 
signatures in the cosmic microwave background (CMB) 
angular distributions. 

All of the above features are generic in models with 
low-scale breaking of neutrino global symmetries, and 
can be understood without reference to a specific model. 
Let us outline the basic physical arguments. 

The traditional seesaw framework relies on the gauge 
quantum numbers of the neutrinos to explain their mass 
spectrum, and leads to a picture where only active neutri- 
nos are light. Global symmetries, on the other hand, may 
involve sterile as well as active neutrinos, and they may 
forbid both Dirac and Majorana mass terms. Neutrino 
masses then appear as a result of spontaneous breaking 
of these symmetries, so that, in this scenario, it is quite 
natural to expect light sterile states. 

If the symmetry breaking phase transition occurs af- 
ter the BBN epoch, both active and sterile neutrinos 
are massless before and during nucleosynthesis. In this 
case, the oscillations, that typically lead to thermal abun- 
dances for the sterile states in the traditional scenario, 
are absent. During BBN, the energy density of the ster- 
ile neutrinos (and of the scalars required to break the 
global symmetries) is determined by their temperature. 
As we show below, this temperature can be significantly 
lower than that of the rest of the cosmic fluid, provided 
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that the rates of certain reactions are sufficiently low. 
This allows our models to evade the BBN constraints on 
AN U . 

Remarkably, the limit from large scale structure on 
the sum of the neutrino masses is also easily avoided by 
the late time neutrino mass models. The breaking of 
the global symmetries gives rise to a set of Goldstone 
bosons, which are coupled to both active and sterile neu- 
trinos. This coupling is sufficiently strong for the sterile 
neutrinos to disappear after they become non-relativistic, 
for example by decaying into an active neutrino and a 
Goldstone boson. As a result, the relic abundance of 
the sterile neutrinos is low, and they do not significantly 
contribute to dark matter despite their large mass. 

Specific Models — There is a wide variety of late time 
neutrino mass models: the neutrinos may be either 
Dirac or Majorana, the number of sterile neutrinos may 
vary, and different choices of global symmetries and their 
breaking patterns can be made. Let us present two sim- 
ple supersymmetric models incorporating the LSND neu- 
trinos. We do not consider models with a single sterile 
neutrino, as they are disfavored by oscillation data. For 
concreteness we construct models with three mass eigen- 
states that are predominantly sterile. We take the global 
symmetry to be U(l) x U(l); a simple possibility that 
allows a heavy neutrino to decay to light neutrino and a 
Goldstone boson. 

Our first theory has three right-handed neutrino super- 
fields, n. There is no overall lepton number symmetry, 
leading to six physical Majorana neutrinos. Above the 
weak scale the theory is described by the superpotential 

W M = W N mssm + W™, 

where Wnmssm is the superpotential of the NMSSM; A, 
A, k and h are coupling constants; and the flavor indices i 
and j run from 1 — > 3. The superfields I, h are the lepton 
and Higgs doublets of the MSSM, s is the electroweak 
singlet field of the NMSSM, and (j>, <j> are the extra elec- 
troweak singlet fields whose vacuum expectation values 
(vevs) give masses to neutrinos. The non-renormalizable 
operators in are generated by integrating out physics 
at scale M; phenomenological constraints discussed be- 
low imply M ~ 10 9 GeV. In theories without an s field, 
the third operator in W^f would be absent, and we would 
expect three light Dirac neutrinos. (If nncj) were allowed, 
the sterile states would be much heavier than the active 
states.) However, for theories such as the NMSSM, where 
the s field acquires a vev of order the electroweak scale, 
the Dirac and Majorana mass terms are of the same order 
of magnitude, vf/M, explaining why the LSND neutri- 
nos are quite close in mass to the active neutrinos. W M 
is the most general superpotential in the neutrino sector 
up to dimension four under the following discrete symme- 
tries: Z3, under which all the fields except <f> and (f> have 



charge 27r/3; Z' 37 under which s, h and h are uncharged, 
q,l,n and </>, <f> have charge 2ir/3, while u c , d c and e c have 
charge — 27r/3; and Z 3 , under which n and <fi both have 
charge 27r/3 while (j> has charge — 27r/3. 

Below the weak scale the renormalizable effective La- 
grangian for the neutrino sector of the theory is 

^ = '/<.;''<"..'•' + gioniUjcj) + h.c. + V(<p, 4>), (2) 

where g = (h) X/M, g = (s) X/M and the scalar potential 
is V = -^ 2 H 2 + K 2 \<f>\ 4 - A 2 H 2 + £ 2 H 4 - (We have 
assumed that SUSY breaking effects generate negative 
soft mass 2 terms for cj), (p.) This theory has two accidental 
U(l) global symmetries: one under which <fi and v are 
charged and another one under which <f>, v and n are 
charged. (While these symmetries are not exact even at 
the renormalizable level, the terms that do not respect 
them are quite small: for example, the cf> 3 term in V((f>) 
is only generated at three loops and can be neglected in 
the present context.) When <fi and <j> acquire vevs, these 
symmetries are broken leading to two pseudo-Goldstone 
bosons G and G, and giving the neutrinos a mass. 

With only minor changes we can construct a theory 
where the six neutrinos are Dirac. There are now three 
singlet left-handed sterile neutrinos, vf, and a total of 6 
right-handed neutrinos n a , coupled via 

W D = Wnmssm +W?, 

W? = X ia kn a h + ^4? +X ia vfn a s ^ + ^ 3 ( 3 ) 

The superpotential W D is the most general up to dimen- 
sion four that is invariant under Z3 x x Z^ (with v s , 
like n, having charges 27r/3 under each Z3) together with 
a lepton number symmetry under which I and v s have 
the same charge and n the opposite charge. 

Below the weak scale the renormalizable effective La- 
grangian for the neutrino sector of this theory is 

C-v = gi a Vin a (i> + gi a f-n a 4> + h.c. + V(4>, 4>). (4) 

The theory has two approximate global symmetries: one 
under which <j> and v are charged, and another under 
which <p an d v s are charged. Again, <f> and (f> vevs 
break these symmetries leading to two pseudo-Goldstone 
bosons and Dirac masses for neutrinos. 

Eqs. (J2J) and imply that g and g are of order 
m u / f, where m v is a scale of order the neutrino masses, 
and / is the scale at which the global symmetries are 
broken. It is important to note that in general the 
couplings of the Goldstones to the neutrinos are not 
diagonal in the neutrino mass basis. Instead, denot- 
ing the mass eigenstates by primes, these couplings are 
of the form {g a fjv' a n'pG + g a pv' a n'pG) (Dirac case) and 

{gapv'aVpG + gapv'av'pG) (Majorana case). 

These theories provide concrete examples of a very rich 
set of theories. A particularly simple theory is obtained 
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by deleting the <f> field and its interactions, and removing 
the Z'z symmetry so that <f) can couple to both doublet 
and singlet neutrino mass operators. In this case there 
is a single flavor diagonal U(l) symmetry and hence a 
single Goldstone, having diagonal couplings to neutrinos 
in the mass basis. 

Constraints — Significant constraints on the parameter 
space of these theories follow from the requirement that 
the total energy density in radiation at the time of BBN 
does not differ significantly from the Standard Model pre- 
diction. This requires that the "hidden sector" fields (</>, 
</>, n and possibly v s , as well as the fermionic partners 
of (j> and <p which will turn out to be quite light) not 
be in thermal equilibrium with the "visible sector" fields 
[y, 7, . . •) before and during the BBN. More precisely, we 
require that the two sectors decouple at a certain tem- 
perature To > 1 GeV, and do not recouple until the tem- 
perature of the visible sector drops below T\y ~ 1 MeV, 
the temperature at which the weak interactions decouple. 
If this is the case, the reheating of the visible sector by 
the decoupling of heavy particles (/i, ir, . . .) and possibly 
by the QCD phase transition will not affect the hidden 
sector. Defining r as the ratio of temperatures of the two 
sectors at the time of BBN, we conclude that the energy 
density in the hidden sector is suppressed by a factor of 
r 4 compared to the naive estimate, and r < 0.3 allows 
one to avoid the BBN constraint even for a very large 
hidden sector. 

The reactions that could recouple the two sectors in- 
clude a 1 <-> 2 process (j> <-> vn, 2 <-> 2 processes such 
as vv <-> nn and vn «-> (jxfi, 2^3 processes such as 
vn «-> 30, etc. Requiring that all these processes be 
"frozen" (r < H) prior to the weak interactions decou- 
pling results in the following constraints on the couplings: 

9ij, 9ia ^ 10~ 5 , gijK, g ta K <; lCT 10 ^ 1 , 

9ij9ij, 9ia9ia £ l(T 10 r- 3 / 2 . (5) 

Note that the coupling k is unconstrained. 

The upper bounds on the coupling g can be trans- 
lated into a lower bound on the symmetry breaking scale 
/. To interpret the LSND result in the model with 
Majorana sterile neutrinos, the low-energy theory must 
possess a mass term of the form my-z^n,, with at least 
some elements of m as large as 0.1 eV. This implies that 
9ijf ~ 0.1 eV, and for a generic flavor structure we obtain 
a bound / J> 10 keV. A similar bound can be obtained 
for the Dirac sterile neutrino case. 

To avoid producing sterile neutrinos by oscillations 
prior to weak interactions decoupling, we require that 
the mass terms mixing active and sterile neutrinos not 
be generated until the temperature of the visible sec- 
tor drops below TV. Scattering in the plasma gen- 
erates "thermal" masses for the bosons, ni 2 (0) ~ 
K 2 n^(T')/T' , where T" is the temperature of the hidden 
sector. The symmetry breaking phase transitions for <j) 



occurs when m 2 (4>) ~ /i 2 . Using /x = /re, we conclude 
that the temperature of the visible sector at the time of 
this phase transition is ~ f/r, implying that f & r MeV 
is necessary for the success of BBN. This in turn imposes 
a lower limit on the couplings, gij > r _1 10 -7 . 

To summarize, BBN considerations lead to a range of 
the allowed values of the scale /, 

10 keV < / < r MeV. (6) 

Considerations of the supernova dynamics may slightly 
raise the lower bound; however, these constraints are 
strongly model dependent 0] . Even though the allowed 
values of / are much lower than the weak scale, the the- 
ory naturally allows for symmetry breaking in this range. 
The only assumption necessary is that </> only feels su- 
persymmetry breaking through its couplings to I and n. 
Then /x 2 is of order <? 2 m| ugY /167r 2 , where msusY is a 
typical soft supersymmetry breaking mass. Since g is of 
order m v j f and / itself is of order /x/k, by eliminating fi 
and g in favor of / and m v we are led to the expression 

, [ml m-susY (r7 . 

For reasonable values of the parameters m„ ~ 0.1 eV, 
"Isusy ~ 100 GeV, k « 10~ 4 , this yields a value of / of 
order 3 MeV, which is quite close to the desired range. 
Analogous considerations apply to the second symmetry 
breaking scale, /. 

Interestingly, the large-scale structure limit on the 
sum of neutrino masses is automatically avoided in 
the models discussed here, and does not lead to addi- 
tional constraints on /. The lower bound on g^ ob- 
tained above implies that the reactions vv <-> nn, (jxf) 
become unfrozen before the sterile neutrinos become non- 
relativistic. These reactions thermalize the hidden sector 
fields with the active neutrinos. The density of thermal, 
sterile neutrinos of mass m s at temperatures T < m s 
is suppressed by a Boltzmann factor e~ ms ^ T ; the excess 
neutrinos disappear either via a decay process n — ► i><p, 
or via an annihilation process nn — * vv. As a result, the 
massive sterile neutrinos do not make a significant con- 
tribution to dark matter. It is only the sum of the masses 
of active, stable neutrinos and the Goldstones that has 
to satisfy the constraints of Ref. . 

Signals in the CMB — The non-standard evolution of 
neutrino masses and densities in our scenario leaves an 
imprint in the CMB inhomogeneities. There are two dis- 
tinct, potentially observable effects Fust, the total 
relativistic energy density at the time of last scatter is 
modified due to the decay and annihilation of the ster- 
ile neutrinos. Second, unlike in the standard cosmology, 
free-streaming of the active neutrinos may be prevented 
by their interactions with the Goldstone bosons. Let us 
consider each of these effects. 

At the time of BBN, the energy density in the hidden 
sector is suppressed. When the reactions vv <-> nn, (jxj), 
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3.78 
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3.70 


3.86 


4.01 


3.91 


4.03 


4.14 
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4.00 


4.11 


4.21 


4.22 


4.29 


4.35 



TABLE I: Effective number of neutrino species during the 
recombination era, N u ,cmb, as determined by the relativistic 
energy density. 



vn «-> <j) become unfrozen, the two sectors thermalize 
and the energy density per degree of freedom in the hid- 
den sector approaches that of the active neutrinos. (The 
active neutrinos themselves are by this time decoupled 
from electrons and photons.) These reactions, however, 
do not change the total relativistic energy density: they 
merely transfer part of the energy from the active neu- 
trinos to the hidden sector states. In contrast, when the 
sterile neutrinos become non-relativistic (T ~ m s ~ 1 
eV) and are depleted by decays and annihilations, the 
total relativistic energy is increased: the depletion pro- 
cess occurs at constant entropy, resulting in an increase in 
temperature as the number of relativistic degrees of free- 
dom decreases. Since the depletion of n's occurs before 
the matter-radiation equality, this will result in a non- 
standard value of the relativistic energy density implied 
by the CMB measurements. In terms of the "effective" 
number of neutrinos JV„,cmb fl^ |. our scenario predicts 



to the Standard Model prediction) is given by 

Al n = 23.3- 13.1 (- g - (3 ~ ns) —] (9) 

where n$ is the number of light neutrinos that are scat- 
tering during the eV era. (It is possible that ns < 3 
if one of the neutrinos is massless or very nearly so, or 
if tog is large enough to make the process z/j — > VjG 
kinematically forbidden for some flavors). Eq. 10 pro- 
vides another experimentally testable prediction of our 
scenario. 

In the theory with no <f> and a single Goldstone, the 
neutrino decays are absent so that scattering can only 
occur via the 2 <-> 2 processes vv <-> GG and vG <-> 
vG. In this case the number of neutrino species which 
scatter is very sensitive to / and to whether the neutrino 
spectrum is hierarchical, inverted or degenerate. 
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, + 2.1hn h jg h 
3 + n G /g u 



Here, g v equals 7/4 for Majorana neutrinos and 7/2 for 
the Dirac case; n^ is the number of the "Higgs" (massive) 
components of the scalar fields responsible for global sym- 
metry breaking that are light enough to be relativistic 
when the reactions vv — > (fxp, vn — > (fxp become unfrozen; 
n s is the number of sterile neutrino species, and no is the 
number of Goldstone modes. (Eq. J5J) includes the contri- 
bution from the superpartners of the symmetry breaking 
scalar fields.) For example, in the explicit models pre- 
sented above, n s — 3 and rig = n^ — 2. Some typical 
values for -/V^cmb are presented in Table [I] For compar- 
ison, while the current sensitivity on N^qmb from the 
WMAP and other CMB analyses is about ±5, the 
sensitivity of the Planck experiment is expected to reach 
the ±0.20 level, providing a test of our predictions. 

Furthermore, at the time of last scatter the mean free 
paths of the light neutrinos and the Goldstones are well 



below the Hubble scale due to the process Vi 



absence of free-streaming leads to a shift in the positions 
of the CMB peaks at large I HHH. This shift (relative 
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